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Abstract — In this paper, we propound a solution named Cog- 
nitive Sub-Small Cell for Sojourners (CSCS) in allusion to a 
broadly representative small cell scenario, where users can be 
categorized into two groups: sojourners and inhabitants. CSCS 
contributes to save energy, enhance the number of concurrently 
supportable users and enshield inhabitants. We consider two 
design issues in CSCS: i) determining the number of transmit an- 
tennas on sub-small cell APs; ii) controlling downlink inter-sub- 
small cell interference. For issue i), we excogitate an algorithm 
helped by the probability distribution of the number of concur- 
rent sojourners. For issue ii), we propose an interference control 
scheme named BDBF: Block Diagonalization (BD) Precoding 
based on uncertain channel state information in conjunction with 
auxiliary optimal Beamformer (BF). In the simulation, we delve 
into the issue: how the factors impact the number of transmit 
antennas on sub-small cell APs. Moreover, we verify a significant 
conclusion: Using BDBF gains more capacity than using optimal 
BF alone within a bearably large radius of uncertainty region. 

Index Terms — Small cells, cognitive radio, Multiuser Multiple- 
Input Multiple-Output (MU-MIMO), Block Diagonalization (BD) 
precoding based on uncertain Channel State Information (CSI), 
auxiliary optimal Beamformer (BF). 

I. Introduction 

Small cells such as femto, pico, and microcells have been 
aroused general interest lately due to their consequential effect 
on enhancing network capacity, stretching service coverage 
and cutting back network energy consumption [1], [2]. 

In many residential and enterprise small cells, according to 
the duration of stay, users can be grouped into two categories: 
sojourners and inhabitants. Besides the duration of stay, other 
features can be pinpointed: the status in the network, the Qual- 
ity of Service (QoS) requirements, the geographical location, 
the regularity or predictability of the entry and exit time. 

Motivated example: A concrete instance in real life can 
well elaborate the above concept. An office area is covered 
by a small cell. The staffs and visitors can be classified as 
inhabitants and sojourners, respectively. Staffs have fixed entry 
and exit time, while visitors enter and depart randomly. Staffs 
have higher status, settled and specific QoS requirements, 
while visitors do not. The staffs and visitors have their own 
area to stay in, such as the working area for staffs and the 
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reception area for visitors. They also have the common area 
to stay like the meeting room. 

Existing approaches and their limitations: Multiple-Input 
Multiple-Output (MIMO) technologies are essential compo- 
nents in 3GPP Long Term Evolution (LTE)-Advanced [3], 
[4]. Among them, multiuser MIMO technology is markedly 
advantageous to enhance the cell capacity especially in highly 
spatial correlated channels [5], [6], [7], since multiuser MIMO 
technology can serve multiple users simultaneously on the 
same frequency. Dirty paper coding (DPC) [8], [9] and Block 
Diagonalization (BD) [10] are two widely adopted techniques 
for multiuser MIMO. DPC is contrived to suit the capacity 
optimal requirement [11], however, it is hard to implement ow- 
ing to its complexity. BD is a practical approach which adopts 
precoding to completely eliminate inter-user interference. The 
deficiency of this technique is that BD precoding and full rank 
transmission requires that the number of transmit antennas is 
not less than the total number of receive antennas. Therefore, 
the number of concurrently supportable users is restricted by 
the number of transmit antennas. 

Cognitive Sub-Small Cell for Sojourners (CSCS) so- 
lution: The above example and the limitations of existing 
techniques inspire us to put forward a solution named CSCS 
in allusion to the following small cell scenario, i) Users can 
be categorized into two groups: sojourners and inhabitants, 
ii) Inhabitants have fixed entry and exit time, while visitors 
enter and depart randomly, iii) Inhabitants have higher status, 
settled and specific QoS requirements, iv) The small cell 
coverage area is relatively large, v) The number of users is 
relatively large, vi) The overlapping area between the inhabi- 
tants' activity area and the sojourners' activity area is relatively 
small. CSCS divides one small cell into two customized sub- 
small cells for inhabitants and sojourners respectively for the 
following reasons. 

1) CSCS is a design in line with the trend of green 
communications [1], [12], which plays a significant role of 
saving energy: 

> Each of two sub-small cells is served by its own sub- 
small cell Access Point (AP). In the case of a small 
overlap between the two groups' activity areas, utilizing 
two sub-small cell APs instead of one single small cell AP 
shortens the average distance between User Equipment 
(UE) and AP. Accordingly, it reduces the total transmit 
power. 

• Considering that inhabitants and sojourners may have two 
different active timetables, CSCS can turn the sub-small 
cell AP off outside the active period of the group which 



it serves. However, AP of a collective small cell can only 
turn itself off when both groups are inactive. 

2) CSCS enhances the number of concurrently supportable 
users. We consider this in the context of BD precoding and 
full rank transmission. The number of transmit antennas on 
AP restricts the number of concurrently supportable users. The 
number of transmit antennas on a single AP is limited by many 
factors, such as the processing speed of AP and the size of AP. 
By utilizing two different sub-small cell APs, CSCS enables 
more transmit antennas to be installed in the whole small cell, 
thereby enhancing the number of concurrently supportable 
users. 

3) CSCS enshields inhabitants by according the priority. 
The QoS requirements and higher status for the inhabitants are 
ensured by a primary sub-small cell. The sojourners are served 
by a secondary sub-small cell. The primary sub-small cell 
and the secondary sub-small cell utilize the same spectrum on 
condition that the level of interference caused by the secondary 
sub-small cell to the primary sub-small cell is kept tolerable 
[13]. 

Main contributions: 

1) Algorithm for determining the number of transmit 
antennas: We propose an algorithm for determining the 
number of transmit antennas on sub-small cell APs helped 
by the probability distribution of the number of concurrent 
sojourners. 

2) BDBF: We proposed an interference control scheme 
named BDBF for secondary systems which are hampered from 
perfect secondary-to-primary Channel State Information (CSI). 
We prove and verify that BDBF can gain more capacity than 
the interference controller using optimal Beamformer (BF) 
alone within a bearably large radius of uncertainty region. 

3) Performance Analysis: i) From simulation and numerical 
results, we find how the factors, such as the standard deviation 
of the duration of stay, the mean of the duration of stay, the 
distribution of arrival rate, and the total number of sojourners 
during the observation time interval, influence the number 
of concurrent sojourners and further impact the number of 
transmit antennas on sub-small cell APs. ii) We show the 
benefit of BDBF. 

The rest of the paper is organized as follows. We present 
the system model in Section II. Section III introduces the 
algorithm for determining the number of transmit antennas on 
sub-small cell APs. Section IV introduces the BDBF scheme. 
We present the simulation results and provide insights on them 
in Section V. Finally, we conclude the paper in Section VI. 

II. System Model 

We consider downlink transmission in a small cell served 
by a single AP. The AP is equipped with Nt transmit antennas 
and every UE is equipped with Nr receive antennas (We listed 
some notations used in the paper in Table I.). The AP adopts 
iVc-subcarrier Orthogonal Frequency Division Multiplexing 
(OFDM), BD precoding and full rank transmission. Based 
on BD precoding and full rank transmission, the number of 
concurrently supportable users is determined by ''Nt/Nr' 1 
[10], [14]. When the number of users in the small cell is greater 
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TABLE I: Notations 



Symbol 


Definition 


r » _l 


the ceiling operator 


Nt 


the number of transmit antennas on the AP 


Nr 


the number of receive antennas at each UE 


AT 


the total number of inhabitants 


AT f+\ 
JV2 (t) 


the total number of sojourners at time t 




the number of inhabitants who stay in the overlap of 




ISSC and SSSC 


the number of sojourners who stay in the overlap of 




ISSC and SSSC at time t 


Q 


the minimum QoS requirements of sojourners 


Q' 


the QoS requirements of inhabitants 


Ntt 


trip iinmhpr nf pnnpiirrpntlv ^nnnnrtarilp sniniimprs 

lllk_ UlllllUd VJ1 k_Ullk_U.ll dill y AUpUUl laLric AUJUUlll^ia 


N st 


the possible number of transmit antennas on SAP 


/V* 

1\ ST 


tVip <jp1pr l l"pH nimhpr o"F tr^nt:mit 'iiitpnn'jt; on QAP 
U1C M_1CI_ICU HUJllL'Cl Ul LliUlMllll alllclllla?! Ull Ortr 


N it 


flip* nntnhpr o"F tr^n<:mit ^ntpnn'K on IAP 
111C 11U111UC1 Ul Ll illlMlllL allLCllllus Ull lr\i 


' f 1 

u ' = {" x) 


the set of inhabitants who stay in the overlap of 




ISSC and SSSC 




the set of sojourners who stay in the overlap of 




ISSC and SSSC at time t 


Us 


the set of sojourners 


Ijv 


the identity matrix of size N X N 


Tr(.) 


the trace operator 


(.)« 


Hermitian transpose 


vcc (X) 


the vector is made up of the columns of matrix X 


® 


the Kronecker product 



than the number of concurrently supportable users, some kind 
of user selection algorithm [14], [15] is implemented. 
We assume for simplicity that 

• Ni inhabitants always stay in the small cell. 

• N 2 (t) sojourners stay in the small cell at time t. 

• Sojourners arrive the small cell according to a time-varing 
Poisson process with arrival rate A (t) and N 2 (t = 0) = 
0. 

• The inhabitants' activity area A\ is a disk of radius R\. 
The sojourners' activity area A 2 is a disk of radius R 2 . 
The small cell coverage area A3 is a disk of radius R3. 
A\ and A 2 partially overlap. A3 accommodate A\ and 

A 2 With R 3 = Rl cos Wl +R 2 cos W2 +R 1+ B^ where ^ ^ 

and the positional relationship of A\, A 2 , and A3 are 
illustrated in Figure 1. 

• The area of Ai, A 2 , and A 3 are, respectively, 5^, Sa 2 , 
and Sa 3 - The area of overlap between A x and A 2 is 

, 5 J 4inA2- 

• Sa 3 is relatively large. S^nAa is relatively small. 

• N 1 inhabitants always stay in the overlapping area be- 
tween A\ and A 2 . 

• N 2 (t) sojourners stay in the overlapping area between A\ 
and A 2 at time t. 

• Inhabitants and sojourners are uniformly distributed in 
their respective activity areas. The number of users is 
relatively large. 

• The inhabitants have higher status. 

A. Cognitive Sub-Small Cell for Sojourners 

CSCS turns the single small cell into two sub-small cells: 
Inhabitant Sub-Small Cell (ISSC) and Sojourner Sub-Small 
Cell (SSSC), each of which is served by its own AP. ISSC 
serves the inhabitants by Inhabitant Access Point (IAP). SSSC 
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Fig. 1 : The system model for CSCS 

serves the sojourners by Sojourner Access Point (SAP). Both 
IAP and SAP exploit BD precoding and full rank transmission. 
The sub-small cells use the same number of subcarriers and 
the same spectrum as that are used by the small cell AP. 
We assume that both sub-small cell APs take inhabitants and 
sojourners in their respective coverage area into account to 
design BD precoding. 

B. Number of Concurrent Sojourners 

Our analysis for SSSC is based on the following assump- 
tions: 

• Sojourners arrive SSSC according to a time-varing Pois- 

son process with arrival rate A (t) and N 2 (t = 0) = 0. 
> The duration of stay of one sojourner is modeled as a 
duration of a continuous time Markov process with m+1 
states from its starting until its ending in the absorbing 
state, where m > 1. The probability of this process 
starting in the state i is a^, where i = 1, . . . ,m + 1. 
The state m + 1 is the absorbing state and the others are 
transient states. We assume that this process will never 
start in the absorbing state, i.e. a m+ i = 0. 
The duration of stay of one sojourner has a phase-type 
distribution [16]. We denote by R the matrix which contains 
the transition rates among the transient states. R is a matrix of 
size m x m. The duration of stay of one sojourner is denoted 
by X. The Cumulative Distribution Function (CDF) of X is 
given by 

F x (x) = 1 - a exp (Rx) 1, x > 0, (1) 

where a = [ati, . . . , a m ], exp (•) is the matrix exponential 
and 1 is an m x 1 vector with all elements equal to 1. 

Remark 1. Under the above assumptions, SSSC can be mod- 
eled as an M t /Gj 00 queue with service time follows a 
phase type distribution, arrival rate A (i) and N2 (t = 0) = 0. 



Authors in [17] have given the probability distribution of 
the number of concurrent users in such queueing system. 
Accordingly, N 2 (t) has a Poisson distribution with parameter 

v (t) = [ A (r) [1 — F (t — t)] dr, t > 0. (2) 
Jo 

The Probability Mass Function (PMF) of N 2 (t) is 

f N2(t) (n) = P (N 2 (t) =n)= exp (-t, (i)) , (3) 

where n G Zq . 

The CDF of N 2 (t) is 

n 

F N 2 (t) (ri) = fN 2 (t) 0) . ( 4 ) 

where n' £ 1^. 

C. Downlink inter-sub-small cell interference 

Due to BD precoding, ISSC and SSSC will not cause 
downlink inter-sub-small cell interference when perfect IAP- 
to-sojourner CSI and SAP-to-inhabitant CSI are acquired by 
IAP and SAP respectively. 

Remark 2. In practice, the perfect acquisition of SAP-to- 
inhabitant CSI is hampered by the lack of explicit coordination 
and full cooperation between ISSC and SSSC, since inhabi- 
tants are scant of willingness to feedback SAP-to-inhabitant 
CSI by using their own bandwidth and power. SAP has to turn 
to blind CSI estimate or other inexact CSI estimates which will 
furnish imperfect SAP-to-inhabitant CSI [18], [19]. Therefore, 
the root cause of interference inflicted by SSSC on ISSC is 
that BD precoding utilizes uncertain CSI. 

In this paper, we consider that IAP-to-sojourner CSI is 
perfect. We suppose the coverage areas of ISSC and SSSC are 
respectively consistent with the inhabitants' activity area A\ 
and the sojourners' activity area A 2 . It is also assumed that the 
inter-sub-small cell interference only exists in the overlapping 
area between A\ and A 2 . 

Since the IAP-to-inhabitant CSI and SAP-to-sojourner CSI 
are perspicuous for IAP and SAP, there is no interference 
inside ISSC or SSSC in our system. In some relevant papers 
about underlay multicarrier cognitive radio systems with un- 
certain CSI, the subcarrier scheduling is the way for steering 
clear of the interference inside the secondary systems [20], 
[21]. 

III. Determining Number of Transmit Antennas on 
Sub-small Cell Access Points 

BD precoding and full rank transmission requires that the 
number of transmit antennas is not less than the total number 
of receive antennas. Therefore, the number of concurrently 
supportable users is limited by the number of transmit anten- 
nas. When the number of users is greater than the number of 
concurrently supportable users, AP will select users to simul- 
taneously serve using some kind of user selection algorithm. 
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QoS requirements and the number of users who simultane- 
ously present make demands on the number of concurrently 
supportable users. 

In our system, we express the required number of concur- 
rently supportable sojourners as N2 (t) Q, where Q indicates 
the minimum QoS requirements of sojourners, Q < 1, and 
N 2 {t)Qe-L+. 

In this paper, we investigate the issue of determining the 
number of transmit antennas under the premise that AP 
can support the installation of selected number of transmit 
antennas. 

We put forward a criterion to determine an appropriate 
number of transmit antennas on Sub-small Cell Access Points: 

• It pursues the minimum number of transmit antennas 
which ensures BD and full rank transmission to be 
validated. 

• It ensures that the mean value of the probability of 
the number of concurrently supportable sojourners being 
adequate during an applicable period is no less than the 
predetermined probability or the growth of this mean 
value of probability led by adding one more supportable 
sojourner is less than a threshold. 

We denote by Njj the number of concurrently supportable 
sojourners. For our system, the probability of the number of 
concurrently supportable sojourners being adequate at time t 



is 



(5) 



P(Nu, t) = F Na{t) (Nu/Q). 

We define the mean value of the probability of the number 
of concurrently supportable sojourners being adequate during 
the time interval [0, , T] as 

E[P(Nu,t)} = ^ [ P(Nu,t)dt. (6) 
1 Jt=o 

According to the first item of our criterion, we have 

N u = N S t/N r Q-N[, (7) 

where Nst is the possible number of transmit antennas on 
SAP and Nu 

Proposition 1. The number of transmit antennas on sub-small 
cell APs can be determined as follows: 

Since Nr and N l are supposed to be known in our system 
model, the selected number of transmit antennas on SAP Ng T 
is determined by (8) or (9): 



N* ST = min{E[P(Nu, t)]> V } 

J\ST 



(8) 



N* ST = min {E [P (N v + 1, t)} -E[P (N Ut t)} < 7} , (9) 

Nst 

where r\ is a predetermined probability and 7 is a threshold. 

We denote by Nit the selected number of transmit antennas 
on IAP. Nit is determined by 



Nit = N R Q' JVi + 



Nl I SA 1 nA 2 



Sa 2 



(10) 



where Njj — N ST /N R Q - N l and Q' indicates the QoS 
requirements of inhabitants. 



IV. Downlink Inter-sub-small Interference 
Control 

As discussed above, for our system, due to the use of BD, 
no inter-sub-small cell interference is caused by downlink 
transmission under the situation that IAP and SAP acquire 
perfect IAP-to-sojourner CSI and SAP-to-inhabitant CSI re- 
spectively. However, the interference inflicted by SSSC on 
ISSC has its origins in BD precoding exploiting CSI with 
inevitable uncertainty. We proposed an interference control 
scheme named BDBF for secondary systems which are ham- 
pered from perfect secondary-to-primary CSI. BDBF utilizes 
BD Precoding based on uncertain CSI in conjunction with 
auxiliary optimal BF. In view of the adoption of OFDM, 
we turn the design of auxiliary optimal BF of BDBF into a 
multi-user multi-subcarrier optimal BF design under channel 
uncertainties. For such optimization problems that contain the 
uncertainty region constraint, S -Procedure is an efficacious 
tool to transform semi-infinite programs into equivalent semi- 
definite programs [19], [20]. We prove and verify that BDBF 
performs better than the interference controller using optimal 
BF alone for gaining more capacity within a bearably large 
radius of uncertainty region. 

A. Block Diagonalization 

For each user u i7 BD algorithm [10] constructs a pre- 
coding matrix W ti „ iiS which achieves the zero-interference 
constraint, i.e. 



H tiU;jiS W tiUiiS = 0, Vj ^ i, 



(11) 



where H t Mj , s G N r x N t denotes the channel matrix from 
transmitter t to user Uj at subcarrier s, N r denotes the number 
of receive antennas at every UE and Nt denotes the number 
of transmit antennas at transmitter t. 

Suppose that N u users are within the coverage area of 
transmitter t and we define 

1 T 



H 



t,Ui,s — 



11 7 



M / 

n t,u i+ i, 



H 7 



Hereupon, Equation (11) has another equivalent form 



H-t,Ui,sWt,Ui,s — 0. 



(12) 



(13) 



From Equation (13), it is deduced that W tjMi;S is a basis set in 
the null space of H ti „ i;S . The Singular Value Decomposition 
(SVD) of H ti „ i;S is expressed as 



*-t,Ui : S 



(1) 



V 



(0) 



t,Ui,S t,Ui,S 



H 



(14) 



where V^l^ e C NtX ( Nt - rank ( ft '.«^)) contains the right- 
singular vectors corresponding to the zero singular values of 
, i.e. Vj°2. s forms a null space basis of H t; „ iiS . There- 



fore W ijBj>s can be constructed by any linear combinations of 
the columns in Vj°t). s . Correspondingly, H t!tli;S Vj°2. s forms 
the equivalent channel. 

The SVD precoding can be performed using the equivalent 
channel: 
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H 



v(0) 

t,ui,s v t,m,s 



U 



t,Ui,S 












V (D y(0) 



H 



(15) 

Where V«, s G C (N t — n k (H ( ,„„ s ))x r an k (H t , Ui , s V<°) ,.) ^ 

tains the right-singular vectors corresponding to the non-zero 
singular values of H tiUi)S Vt°2 s . a is used as SVD 

precoder and U^ u s is used as SVD decoder. Consequently, 
the precoding matrix W ()Uji8 can be chosen as V^ -}. gV^. g . 
For the full rank transmission, i.e. rank (jlt tUiiS 'V[ ^. = 

min |iV r , , _/V t — rani; ^Hj iUijS ^ |, in order to guarantee that 

the null space of H t; „ iiS is not empty, the condition that N t > 
N r N u is required. In this paper, we consider the situation that 
the number of streams sent to a UE equals the number of 
receive antennas, i.e N t = N r N u . Accordingly, W( iUji3 is a 
matrix of size N t x N r . 



B. BDBF 

As explained in Remark 1, SAP-to-inhabitant CSI is un- 
able to be perfect. The SAP-to-inhabitant channel matrix at 
subcarrier s can be expressed as 



(16) 



where HsAPji.s is the estimated channel matrix from trans- 
mitter SAP to user J, G Uj at subcarrier s and AHs^^/.j is 
the channel uncertainty matrix from transmitter SAP to user 
Ii G Uj at subcarrier s. In this paper, we adopt a channel 
uncertainty model [22] which defines the uncertainty region 



as 



A(e) = {AH SAPJiiS Tr{E}<e 2 } 



,AH" 



(17) 



where E = AH SA pj itS L^i± SAP I . s . 

Based on the assumption that the inter-sub-small cell inter- 
ference only exists in the overlapping area between A\ and A 2 , 
SAP only need to take the inhabitants in the overlap between 
Ai and A 2 into account for restraining interference. 

Proposition 2. For secondary systems which are hampered 
from perfect secondary-to-primary CSI, secondary systems 
can exploit BDBF: BD precoding based on uncertain CSI in 
conjunction with auxiliary optimal BF to control interference 
to primary systems. 

After BDBF, the transmitted signal from SAP to sojourner 
Si can be expressed as 

SsAP,S t ,s = W SAP,S t ,sP SAP,S t ,s*-SAP,Si,s (18) 

where 

X-SAP,ii,s G C NstX1 is the information symbol vector trans- 
mitted from SAP to Ii at subcarrier s with covariance matrix 
^ {^-SAP,i i , s x'sAP,i i , a } = Ijv| t , 

Psap,i z ,s G C NrXN st is the auxiliary optimal BF, 

and WsAP,Si,s is the BD precoding matrix of size Ng T x Nr. 



WsAP,Si,s is generated 



H 



SAP, Si, a 



tt7 

n SAP,Ji,s ■ 



tt7' 

• n SAP,I 



by 

T 



H 



SAP,Si,s 



ttV tt7 tt7' 

"-SAP.Sls ■ ■ ■ n SAP,Si_i,s n SAP,S i+1 ,s ■ 



and H SAPjS . jS 



tt7' 

■ t± SAP,S N2{t) , 



For inhabitant f G U I , the received signal at subcarrier s 
is given by 



— Uj^ Pi 7 . gH/AP./j , s W IAP,Ii,s x IAP,Ii,s 



(19) 



where 

*-iAP,ii,s G C JV/TXl is the information symbol vector trans- 
mitted from IAP to Ii at subcarrier s with covariance matrix 

E {'X-IAP,I i ,s^IAP,I i ,s} = IiV/r' 

n/ ijS G C^"* 1 is the i.i.d. circularly symmetric white Gaus- 
sian noise at subcarrier s at reciever /; with covariance matrix 
E{n It , a n¥ iyS } = I Nr , 

~WiAP,i i;S is the BD precoding matrix of size Nit x Nr, 

and UYap ,i f , s ls tne SVD decoding matrix of size Nr x Nr. 

W/ap,;„ s and Uf AP j. s are generated by H/ap,/„ s = 



H7 



IAP,I i}S H JAP>iSliS . . 



HT 



' JV2O 

HT 



and H/AP,ij, a 



Theorem 1. 77ze auxiliary optimal BF of BDBF can be formu- 
lated as the following multi-user multi-subcarrier optimal BF 
design under channel uncertainties: trying for the maximum 
capacity of secondary system under the precondition that the 
interference inflicted to the primary system is kept below the 
threshold. 



(PI) max £ £ logdct 1^ + -— A (20a) 

QsAP, Si ,s s=1 s , eUs \ iV ST J 

s-t. QsAP, Si ,s > 0, VS t G U S , Vs, (20b) 



SitUs 



Tr (QsAP, Si , s ) < Psap, VS 4 G U S , Vs, 



(20c) 



Tr(B) < C, VAH S , PJ „ S G A (e) , VS 4 G C/ s , VJ* G U T , V«, 

(20d) 

QsAP.Si.s = E {PsAP,S„sPsAP,Si, s } ' 
QsAP.Si.s = WsAP.Sj^QsAP.Si.sWMP.S,,^ 
A = A.SAP,S i ,aQsAP,S i ,s^SAP,S i ,sj 

B = HsAP,f„ s Q Sy ip,s i , s HsAP,;„ S i 

Psap is the maximum downlink transmit power of SAP, 
and £ is the inter-sub-small interference tolerable threshold 
on each subcarrier for inhabitant Ii G Uj. 

Q SAP s- s i s a positive semidefinite matrix to render certain 
that the elements on the main diagonal of A and B are 
real and non-negative. Additionally, QsAPSi s ™ a symmetric 
matrix to ensure that linear matrix inequality (LMI) holds. 
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A (e) is a set of infinite cardinality and accordingly the 
inequality (20d) imposes infinite constraints. This means that 
(PI) is a semi-infinite program. (PI) can be transformed into 
an equivalent semi-definite program (P2) with the constraint 
(2 Id) which contains finite constraints instead of the constraint 
(20d): 



H? 



N c 



(P2) o max ]T logdet(l w „ + ^-A) (21a) 



SNR 



" s—lbiGUs 

s-t. Qsap^s > 0, y,S l G U S , Vs, (21b) 
E Tr (0sap,s () ,) < p sap, VSi G C/ s , Vs, (21c) 

Si&Js 



aI NTXNa + G J n 

J k - ae 2 

Ma > 0, Wj G t7j, VS* 4 G U s , Vs 



> 0, 



(21d) 



Proof: According to the relationship be- 
tween the trace operator and the vec operator 
Tr(C w DC) = 



vec {C) n (I <g> D) vec (C) and 



Tr (C«D 



vec (C) n vec (D), the inequality (20d) 



can be converted to an equivalent form: 



-vec (ah|mp >JiiS ) (l NR <g> Qsap,s z ,s) vec (aH^,/.,,) 

-2% (vec (Q SAPtSii3 H il% APJi ^ n vec (AH^ P ,, j:S )j 

-Tr ^Hs^p,i i , 3 Os J iP,s i , s Hs / ,p i ,. i ,j 

+C >0, 

V/, G t/j, VSi G J7s, Vs (22) 

The result of S -Procedure [23], [24] presents a condition 
that makes the infinite constraints imposed by the inequality 
(22) hold: if and only if there exists a ^ such that the 
inequality (23) is ture. 



C&NtxNr + G 



J* 

k — ae 



> 0, V/i G [/,, VS, G (7 S , Vs 
(23) 

In r ® QsAP,Si,s> J 

fc = C - 



where G = 

-vec(Qs AP;Si;S H^ 4 p ;7 . s ), and 

Tr (HsAP.Ji.sOsAP.Si.sHsAP./i,*)- ■ 

The semi-definite program (P2) can be solved by interior 
point methods [24], [25]. 

Theorem 2. BDBF can gain more capacity than the interfer- 
ence controller that using optimal BF alone within a bearably 
large radius of uncertainty region. 

Proof: We denote by W a SAPSt S and W b SAPSt S 
the BD precoding matrix which is generated by 



H 



SAP, Si, s H Jj4 p JliS . . 



.hi;, 



1 . 



and 



the BD precoding matrix which is generated by 

h sap,s,, s = [h<?ap,s,, s ] , respectively. Since W SAP<S . >a is 
a basis set in the null space of tlsAP s- s> tnen f° r ^ — L ' 
we have Tr (HsAPj i ,.Vr%AP,s i ,.W&p,s i ,.tisAP,i i ,.) < 
Tr(H M P, /j , 3 W| AP , s . s W^p, s . s H^ Pi/ . s ), where i 
is a constant. Accordingly, for e 2 < l, when Tr(B a ) = 
Tr (B fc ) = C, we have Tr (Q SA P, Si , s ) > Tr (Q^ P , Sj j, where 
B" = Hmp./^W 



Tb 



sQsAP,S i ,s^SAP,S i ,s'Hs 



H AP,h,s and 



B 6 - HsAP,I i ,a^V b SAP,S i ,aQ b SAP,Si,s^^AP,S i ,s' H -'sAP,I i ,s- 

Therefore, for e 2 < i, based upon Theorem 1., the capacity 
of SSSC using W% APS . s is larger than that of SSSC using 

W SAP,S„s- ■ 

We simulate the SSSC using ~Wg APS s and Wg APSj8 for 
different e 2 in Section IV. We find that i can be much larger 
than the normally acceptable radius of uncertainty region. 



V. Simulation 

A. Number of Transmit Antennas on Sub-small Cell Access 
Points 

Since the number of transmit antennas on SAP and that 
on IAP are in direct relation with the number of concurrent 
sojourners, we design 8 data sets to observe and analyze the 
number of concurrent sojourners in 8 representative circum- 
stances. 

We consult a medium-sized software enterprise in Paris 
to help us to design data sets to approximate the actual 
situation. We receive the following information on visitors in 
this enterprise: 

• There are more than fifty visitors daily at an average. 

• In view of long term observation, there is no such a 
working day when the number of visitors is particularly 
large or small within five working days. 

• After the holidays, the number of visitors is inclined to 
peak. 

First, we exclude those days when the number of visitors 
is inclined to peak after holidays. Based on the second item 
above, we consider that the observations of the number of 
arrivals at the reference time on different working days are 
independent and identically distributed (i.i.d.). We observe the 
sample mean of the cumulative probability of the number of 
concurrent sojourners during one entire working day based 
on the first six data sets. Here, we set Q = 1. We adopt 
E[P(Nu, t)] which is defined in Equation (6) to indicate 
the sample mean of the cumulative probability of the number 
of concurrent sojourners during the observation time interval. 
We interpret P [Ny , t) as the probability of the number of 
concurrent sojourners being less than Nu at time t. 

Then we consider separately those days after holidays. We 
also consider that the observations of the number of arrivals 
at the reference time on different working days during those 
days are i.i.d.. We observe the sample mean of the cumulative 
probability of the number of concurrent sojourners during one 
entire working day based on the last two data sets. 

We list the attributes of these 8 test data sets in Table (II). 
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TABLE II: Attributes of test data sets 



Data set index 



Observation time interval (9 am to 5 pm) 



Total number of sojourners during the observation time interval 



Sample mean of the duration of stay of one sojourner (minute) 



Standard deviation of the duration of stay of one sojourner 



Arrival rate (arrivals per time slot) 



Data set 1 


Data set 2 


Data set 3 


Data set 4 










60 


60 


60 


60 


60 


60 


90 


90 


3.7947 


148.0447 


3.7947 


148.0447 


60/48 


60/48 


60/48 


60/48 


Data set 5 


Data set 6 


Data set 7 


Data set 8 










60 


60 


90 


90 


60 


60 


60 


60 


3.7947 


148.0447 


2.7809 


134.5144 


5(10-11 am), 30/42(else) 


5(10-11 am), 30/42(else) 


60/48 


60/48 



The staffs' working hours are from 9 am to 5 pm. The vis- 
itors arrive randomly during working hours. The observation 
time interval is set to be consistent with the working hours 
for all data sets. We divide the observation time interval into 
48 time slots, each 10 minutes in length. Data set 1 and data 
set 2 have the same attributes except the standard deviation 
of the duration of stay. These two data sets are considered as 
the benchmark for comparison and designing the other data 
sets. Data set 3 and data set 4 are designed by adding 30 
minutes to every sample of the duration of stay in data set 1 
and set 2. Data set 5 and data set 6 are designed by changing 
the immutable arrival rate of data set 1 and data set 2 to 5 
arrivals per time slot during 10-11 am and 30/42 arrivals per 
time slot during the other working hours. Date set 7 and set 8 
are designed by adding 30 to the total number of sojourners 
during the observation time interval of data set 1 and data set 
2 and make the the standard deviation of the duration of stay 
of data set 7 and data set 8 to be close to that of data set 1 
and data set 2 respectively. 

We exploit EMPHT algorithm [26], [27] to fit the phase type 
distribution to the data of duration of stay in each data set. 
The EMPHT is an iterative maximum likelihood estimation 
algorithm which outputs the m X 1 vector a and m x m 
matrix R. We chose the number of phases as 4, i.e. m = 4. 
We acquire these two parameters based on each data set after 
sufficient iterations when the likelihood function reaches zero 
growth and accordingly obtain the CDF of the duration of stay 
for each data set. 

We further calculate the mean value of the cumulative 
probability of the number of concurrent sojourners during the 
observation time interval base on each of these 8 data sets 
and showed them in Figure 2. From them, we can observe 
how the following factors influence the number of concurrent 
sojourners and further impact the number of transmit antennas 
on sub-small cell APs: 

• the standard deviation of the duration of stay 

• the mean of the duration of stay 

• the distribution of arrival rate 

• the total number of sojourners during the observation time 
interval 



i 

0.98 
(1.96 
0.94 

5 0.92 

6 0-9 
^ 0.88 

0.86 
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(d) 

Fig. 2: Sample mean of the CDF of the number of concurrent 
sojourners during the observation time interval base on each 
data set 



x 



From Figure 2, we observe that as Njj increases, its sample 
mean of the cumulative probability of the number of concur- 
rent sojourners during the observation time interval based on 
8 data sets all exhibit uptrend and the rate of growth declines 
gradually. The values of E [P (Nu = 10, t)] based on data set 
1 until data set 8 are: 0.8237, 0.9354, 0.4711, 0.5223, 0.7550, 
0.8849, 0.4048 and 0.6905. The values of E [P (N v = 20, *)] 
based on data set 1 until data set 8 are: 0.991 1, 0.9719, 0.9841, 
0.9857, 0.8908, 0.9554, 0.9852 and 0.9895. The pair of 
n = win Nu {E[P(Nu + l, t)]-E [P{N V , t)\ < 0.001} 
and E [P (n, t)\ based on data set 1 until data set 8 
are: (19, 0.991), (15, 0.9712), (24, 0.9912), (24, 0.9911), 
(38, 0.9741), (19, 0.9553), (26, 0.9961) and (22, 0.9912). 

From the first two sets of experimental results, we can 
observe that: the number of concurrent sojourners during the 
observation time interval being more than 20 has a very small 
probability (<0.03) for all data set 1, 2, 3, 4, 7 and 8 and a 
small probability (<0.11) for data set 5 and 6. This number 
being no more than 10 has a large probability (>0.8) for data 
set 1 and 2, followed by data set 5 and 6 (>0.75). For data 
set 3, 4, 7 and 8, This number being no more than 10 and 
between 10 and 20 have almost the same probability (« 0.4). 

Referring to Figure 2, we can find that the increasement 
of the mean of the duration of stay, non-uniform arrival rate 
or the increase of the total number of sojourners during the 
observation time interval improves the probability of the num- 
ber of concurrent sojourners being a relatively large number. 
The increasement of standard deviation of the duration of stay 
brings down this probability. However, when the mean of the 
duration of stay is relatively large, the standard deviation of 
the duration of stay has minimal impact on this probability as 
shown in Figure 2 (b). Accordingly, the number of transmit 
antennas on sub-small cell access points are chosen to be 
relatively large when the mean of the duration of stay or the 
total number of sojourners is relatively large, the arrival rate is 
non-uniform or the standard deviation of the duration of stay 
is relatively small. 

From the third set of experimental result, we can see that 
the augmentation of standard deviation of the duration of 
stay causes that the convergence of the sample means of the 
cumulative probability of the number of concurrent sojourners 
during the observation time interval to be more slowly and 
the low growth rate emergences earlier. But when the mean of 
the duration of stay is relatively large, this phenomenon is not 
obvious. So we need to choose a relatively small value for rj in 
Equation (8) or 7 in Equation (9) when the standard deviation 
of the duration of stay is relatively large and the mean of the 
duration of stay is relatively small. 

B. BD preceding based on uncertain CSI in conjunction with 
auxiliary optimal BF for downlink inter- sub- small interference 
control 

We consider that 2 sojourners and 2 inhabitants stay in 
the overlap of ISSC and SSSC, where all the sojourners and 
inhabitants have the same distance from the SAP. The path loss 
component is normalized to one. HsAP,Si.s and HsAP,it,a are 
both modeled as a matrix with coefficients which are i.i.d. 



£.48 







i i ■ - -* J 


* " ~ ~J 






























BD includes inhabitants (perfect CSI) 
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Interference threshold (W) 

Fig. 3: Average capacity of SSSC vs interference threshold 




Interference threshold (W) 

Fig. 4: Average interference on each subcarrier vs. interference 
threshold 



circularly symmetric, complex Gaussian random variables, 
with zero mean and unit variance. AHsap,i,,s is modeled as a 
matrix with coefficients which are i.i.d. circularly symmetric, 
complex Gaussian random variables, with zero mean and 
variance a 2 . We set that e 2 in Equation (17) equals to a 2 . 
The number of receive antennas at every UE is set to be 
2 and the maximum Signal-to-Noise Ratio (SNR) at every 
UE is set to be 20 dB. The maximum downlink transmit 
power of SAP is 10W. The SAP exploits 4 subcarriers for 
downlink transmission. Each subcarrier occupies a normalized 
bandwidth of 1 Hz. 

For the purpose of comparing the performance of BDBF 
with the interference controller that uses optimal BF alone, 
we design two test systems: 1. downlink BD precoding for the 
secondary sub-small cell takes both sojourners and inhabitants 
who stay in the coverage area into consideration, i.e. the SSSC 
using Wj^p s . ; 2. downlink BD precoding for the secondary 
sub-small cell takes only sojourners into consideration, i.e. the 
SSSC using W| AP s . s . We simulate the above systems which 
exploit BD precoding based on SAP-to-inhabitant CSI with 
different radius of uncertainty region (e 2 =0, 0.1, 1, 5 and 10) 
and fed perfect SAP-to-inhabitant CSI to the optimal BF. We 
plot the average capacity of SSSC and the average interference 
to inhabitants on each subcarrier which are obtained by using 
2000 Monte-Carlo realizations in Figure 3 and Figure 4 
respectively. 

Figure 3 reveals the following information. 1. For the 
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system adopts BD precoding including inhabitants, the average 
capacity of SSSC shakes off the shackles of interference 
threshold within the small radius of uncertainty region. For 
e 2 = 0.1, the average capacity of SSSC achieves 54.7651 
bps when the interference threshold is 1.75W. After that it 
oscillates on small scale and no longer increases with the 
enhancement of threshold. This is because BD precoding 
based on uncertain CSI can inhibit the interference below the 
interference threshold within the small radius of uncertainty 
region. Figure 4 exhibits that the average interference to 
inhabitants on each subcarrier for e 2 = 0.1 is maintained 
below the threshold when the threshold is no less than 0.25W. 
For e 2 = 0, since BD precoding based on perfect CSI, the 
interference is completely eliminated and accordingly the aver- 
age capacity of SSSC has no relationship with the interference 
threshold. It fluctuates slightly with the fluctuation of channel 
between SAP and sojourners. For e 2 = 1, 5, and 10, the 
auxiliary optimal BF plays the leading role in inhibiting the 
interference. The actual interference is always equal to the 
threshold. The average capacity of SSSC exhibits uptrend as 
the threshold is enhanced. 2. System adopts BD precoding 
including inhabitants achieves higher capacity compared with 
that exploits BD precoding excluding inhabitants within a 
considerably large radius of uncertainty region. The average 
capacity of SSSC of the first system declines as the radius 
of uncertainty region increases. It should be noted that until 
e 2 increases to 10, the average capacity of SSSC of the first 
system is still higher than that of the second system with 
e 2 = 0. This verifies Theorem 2. 

VI. Conclusion 

In regard to design issues in CSCS, we put forward practical 
propositions for determining the number of transmit antennas 
and controlling downlink inter-sub-small cell interference. We 
found how the factors influence the number of concurrent 
sojourners and further impact the number of transmit antennas 
on sub-small cell APs. We proved and verified that BDBF 
performs better than the interference controller using optimal 
BF alone for gaining more capacity within a bearably large 
radius of uncertainty region. 
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